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Abstract 

The Hamiltonian of persistent current qubit is found within well known quantum mechanical 
procedure. It allows a selfconsistent derivation of the current operator in a two state basis. It is 
shown that the current operator is not diagonal in a flux basis. A non diagonal element comes 
from the finite inductance of the qubit. The results obtained in the paper are important for the 
circuits where two or more flux qubits are coupled inductively. 

PACS numbers: 03.67.Lx, 85.25.Cp, 85.25.Dq, 85.35.Ds 



I. INTRODUCTION 



Josephson-junction qubits are known to be candidates for scalable solid-state quantum 
computing circuits [1] . Here we consider a supercobducting flux qubit which has been first 
proposed in [2] and analyzed in and j^. The qubit consists of three Josephson junctions 
in a loop with very small inductance L, typically in the pH range. This insures effective 
decoupling from the environment. However, in the practical implementation of flux qubit 
circuitry it is important to have the loop inductance as much as possible consistent with 
a proper operation of a qubit. A relative large loop inductance facilitates a qubit control 
biasing schemes and the formation, control and readout of two-qubit quantum gates. These 
considerations stimulated some investigations of the role the loop inductance plays in the 
operation of a flux qubit 0], Q. The main goal of these works was the calculation of 
he corrections to the energy levels due to finite inductance of the loop. In the early work 



5, 



these corrections have been obtained by perturbation expansion of the energy over small 



parameter (3 = L/Lj, where Lj is the Josephson junction inductance. The extension to 
large /3's (up to (5 ~ 10) had been considered in |7|. However, it is important to realize 
that for finite loop inductance the interaction between two state qubit with its own LC 
circuit cannot in general be neglected. If f3 is not small, as in [7J, this interaction can have 
substantial influence on the energy levels. Unfortunately, this interaction in []| has been 
completely neglected. 

In principle, the account for a finite loop inductance (even if it is small) requires the 
correct construction of quantum mechanical Hamiltonian of a qubit, which contains all 
relevant interactions. This has been done in [61, where the effective Hamiltonian has been 

n 

obtained by a rigorous expansion procedure in powers of j3. As was shown in |6[, one of the 
effect of the interaction of a flux qubit with its own LC oscillator is the renormalization of the 
Josephson critical current. The inclusion of circuit inductances in a systematic derivation 
of the Hamiltonian of superconducting circuits has been done in W. It allows the correct 

q n 

calculations of the effects of the finite inductance both for flux [2[ and charge [10J qubits. 

In this paper we investigate another physical effect which comes from finite loop induc- 
tance. Namely, we show that the finite loop inductance results in the additional term of the 
current operator in the flux basis: 

T=At z + Bt x , (1) 
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where Tz and Tx are Pauli matrices in the flux basis. The quantities A and B in (JIJ are 
calculated in the paper the B being conditioned by the finite loop inductance: for L = the 
second term in (JTJ) is absent. Though for the usual qubit design with small loop inductance 
this second term is relatively small, nevertheless, it might give noticeable effects for large /?'s 
for the arrangements when two flux qubit are coupled either via a common inductance ]j or 
inductively coupled via a term MI1I2 in the Hamiltonian, where M is a mutual inductance 
between qubit's loops, Ii, I2 are the current operators of the respective qubits. 

The paper is organized as follows. In Section II we start with the exact Lagrangian of a 
flux qubit with finite loop inductance. In section III with the aid of well known procedure 
we derive rigorously the quantum qubit Hamiltonian. The current operator is studied in 
Section IV, where we show that in general it is not diagonal in the flux basis. The matrix 
elements for the current operator are calculated in Section VI, where in order to obtain 
analytical results we consider a flux qubit with a small loop inductance. 

II. LAGRANGIAN FOR THE FLUX QUBIT 

We consider here a well known design of the flux qubit with three Josephson junctions 
jsL A) which is shown on FigHJ 

Two junctions have equal critical current I c and (effective) capacitance C, while those 
of the third junction are slightly smaller: al c and aC, with 0.5 < a < 1. If the Josephson 
energy Ej = J c $ /2vr is much larger than the Coulomb energy E c = e 2 /2C, the Josephson 
phase is well defined. Near $ x = $o/2, this system has two low-lying quantum states 3,0. 
The Lagrangian of this qubit is the difference between the charge energy in the junction 
capacitors and the sum of Josephson and magnetic energy: 

1=1 1=1 

where Vi is the voltage across the junction capacitance Ci, which is related to the phase 
4>i by the Josephson relation Vi = ($o/27r)0j; $ is the flux trapped in the loop: 

i=l 
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FIG. 1: A flux qubit, where an external magnetic flux <I>x pierces the superconducting loop that 
contains three Josephson junctions and inductance L. Two Josephson junctions are considered to 
be identical, Ej\ = Ej2 = Ej, C\ = C2 = C ', and Ej3 = aEj, C3 = aC. 



Next we make the following definitions: <p = <pi + <p2 + 4>3, 4>i + 4>2 
terms of these new phases Lagrangian O) takes the form: 



20, 



2 X . In 



h 2 



16E C 

where j3 = 2kLI c /§o. 



oc- 



16K 



E, 



-03 + 2Ej cos 9 cos x + o>Ej cos (if - 29) - — (<p - tp x ) (4) 
c 



III. CONSTRUCTION OF HAMILTONIAN 



n <p = T— = a —[V-29) (5) 



Conjugate variables are defined in a standard way: 

IdL _ h 

hdfi ~ a 8E c 

1 dL h ■ h ( „ -\ 

n ° = nTe=^ d - a ^- 2d ) (6) 
idL h . 

From these equations we express phases in terms of conjugate variables: 

. %E c 2a + l 8E C fo . 

* = ~h + (8) 
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(9) 



h 

In terms of conjugate variables Lagrangian (jlj) takes the form 



X=^n x (10) 



L = 2E c n 2 + 2E c n\ + AE c ^-^n 2 + %E c n 9 n v (11) 

a 

E j 2 

+2Ej cos 9 cos x + cos (ip — 29) — (if — <fx) 

Zp 

The Hamiltonian is constructed according to the well known rule: 

H = Hritpip + hrie9 + hn x x — L (12) 

Finally we obtain: 

H = AEc^^nl + 2E c n\ + 2E c n 2 e + 8E c n n v + U { X , 9, ip) (13) 
a 



where 



E 

U( X , 9, if) = -2Ej cos 9 cos X - uEj cos (<f - 29) + -| (<p - <Px? (14) 



Hence, equations of motion for the phases (jBj), (jUJ), (fTU|) are simply = 9 
\§n' : > X = IJrF' ^ ne e Q ua tions of motion for conjugate variables are: 



IdH aEj . Ej 



sm(ip -29)--±(ip- f X ) (15) 



IdH 2aEj . IE j 
n e = --^-QQ= —sm(ip-29) — sin cos x (16) 



n x = --— = -— cos£sin X (17) 

Below we consider Hamiltonian (jTHj) as quantum mechanical with commutator relations 
imposed on its variables 

[<p,n v ]=i; [9,n e ] = i; [ x ,n x ]=i (18) 
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IV. CURRENT OPERATOR 



From the first principles a current in the loop is equal to the first derivative of the state 
energy relative to external flux: 

d$ x 

This expression can be rewritten in terms of exact Hamiltonian of a system: 

/ = („! — 1„) (20) 
From (|2U|) we would make ansatz that the current operator is as follows: 

■ 9 " (21) 



However (|21jl is not a consequence of (|2()jl . Therefore, the ansatz (|21jl must be proved in 
every case, since the current operator in the form of Eq. (J2TJ) has to be consistent with its 
definition in terms of variables of Hamiltonian H. The prove for our case is given below. 

The current operator across every junction is a sum of a supercurrent and a current 
through the capacitor: 

It = I sin <p t + ^-CCpi (2 = 1,2) (22) 

J 3 = al sin yj 3 + a ^ C< Ps ( 23 ) 

Since the current in a loop is unique the equations (J22J) and ()23j) must give identical result. 
This is indeed the case if we express phases = 1,2,3) in terms of <p, 9, % an d use the 
equations p. p. (flftl). (fISl). (fTHjl . (fT7jl . For every I- in (|22^ . (|2^)l we obtain the same 
expression 

/ = "/.^ (24) 

which is independent of parameters of a particular junction in the loop. From the other 
hand the expression (|2*^|) can be obtained from our Hamiltonian (fT3*j) with the aid of (|2ip. 
Therefore, the equation (|21j) gives us the true expression for the current operator. It is 
important to note that the proper expression for the current operator (|24jl cannot be obtained 
without magnetic energy term in the original Lagrangian (J2J). 

It follows from (|T3*|) and (}2"4"j) that J, H ^ 0. Therefore, an eigenstate of H cannot 
possess a definite current value. 



A. Current operator in a two-state basis 



Suppose a system is well described by two low lying states \^±) with corresponding 
eigenenergies E±: 

H\V±) = E±\V±) (25) 
Within this subspace Hamiltonian can be expressed in terms of Pauli matrices ax, &y, &z' 



fr E + + E„ E + -E_ 
tl = z o z 



(26) 
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with a z \^±) = T|*±)- 

Now we calculate the matrix elements of the current operator (J21|) within this subspace. 
According to (|19j) . (j2(J|) and (J21j) diagonal matrix elements are: 

dE± 



(*±|I|*±) 



<9$ 



x 



(27) 



In order to find nondiagonal matrix elements of the current operator we use the expression 



OH 

(n\ -qT- \n') = (E n , - E n )(n 



dn' 



(28) 



which is obtained by differentiating of the identity (n \ H \ n') =0 with respect to parameter 
A. Hence, the nondiagonal elements of the current operator are: 



(*_| I\^+) = (#+| / = (E + - EJ) ( * 



<9$ 



x 



Therefore, we can express the current operator in terms of Pauli matrices: 



d {E + + E_ 



I 



d (E+-E_ 



d$x \ 2 J d<S>x 
where I is the unity matrix. 

Below we consider two low lying states of a flux qubit 



a z + (E+ - E_ 



<9$ 



x 



ox 



(29) 



(30) 



E± = E ± Ve 2 + A 2 



(31) 



where Eq and the tunneling rate A are independent of the external flux $x, and the quantity 
e is linear function of the flux, e = EjXfx, where A is a numerical factor which depends on 
qubit parameters a and g = Ej/E c , fx = &x/&o — 1/2. 
Therefore, for the flux qubit we get in eigenstate basis: 



H = E - A e a z 
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(32) 



dA F 



<9$ 



-a z 



2A £ ( 



x 



3$ 



A 



ox 



(33) 



where A e = ^/e 2 + A 2 . 

Transformation to the flux basis is obtained via the rotation around y axes in a two level 
subspace with the aid of the matrix R = exp (i£oy/2), where cos£ = e/A e , sin£ = A/A £ : 
Rr 1 (7zR = T z cos^+rx sin£, R~ 1 a x R = —t~z sin £+Tx cos £, where tx, are Pauli matrices 
in a flux basis. Hence, we get for Hamiltonian (J32|) and current operator (J33|) in the flux 
basis: 

H = -sr z - At x (34) 



I = - 



dA £ e 
d§x A e 



2A, ( 



Therefore, the current operator is not diagonal neither in the flux basis nor in the eigen- 
state basis. 

The stationary state wave functions can be written as the superpositions of the wave 
functions in the flux basis, ^l^r where L, R stand for the left, right well, respectively: 
\l/± = a±$ l + b±^R, where 

A , eTA e 



a± 



(36) 



2A e (A e =p e) ' ^ v/2A £ (A eT 5) 
The coefficients a±, b± are defined in such a way, that t z \^l) = — t z\^r) = + |^fi)> 
t~x\^l) = +\^r), Tx\^?r) = In terms of the functions ^l,^r the cross term 



will read 
5*, 



6>$ 



96, 



9$ 



A 



9$ 



<9$ 



6_6 , ( # 



9*, 



9$ 



A 



(37) 



The results obtained up till now are exact in that we did not make any approximation to 



the Hamiltonian (|13|). However, in order to calculate e and cross term (\EL 
(JHEJ) we need some approximate procedure. 



in flUD, 



V. APPROXIMATION TO QUANTUM MECHANICAL HAMILTONIAN 



In order to calculate the matrix elements of the current operator we have to find the wave 
functions of two lowest levels of Hamiltonian ([13)1 . First we single out of the potential (JUjl 
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the fast variable ip, which describe the interaction of the qubit with its own LC circuit. The 
point of minimum cpc of U(x, 9, <p) $111 with respect to (p is defined from dU/d<fi = 0: 

<-Pc = ¥x ~ a(3 sin (ip c - 29) (38) 

In the vicinity of <pc the potential U (x, 9, V 5 ) can be written as: 

U ( X , e,cp)*U ( X , 8, <p c ) + ^ 2 + $ 2 ^- cos {ipc - 29) (39) 

where <p is a small operator correction to ipc'- f — fc + <P- 

As is known the potential U(x,9,<po) has a degenerate point at $ x = $o/2. Assuming 
/x « 1, « 1 we obtain near this point: 

<f C = 7T + 2nf x - a(3 sin 29 (40) 

From (|2*4"j) and (j4*Uj) we find a current operator in "coordinate" representation: 

J = / asin2# (41) 

For U (Xi ®i Pc) we obtain near degeneracy point 

o?3Ei 

U (x, 0, <p c ) = ~2Ej cos 9 cos x + ®Ej cos 29 + a2nf x Ej sin 29 sin 2 29 (42) 

Below we follow the procedure described in At f x = 0, the potential (f4*2"J) has two 

minima at x = 0, 9 = ±9*, with cos#* = l/2a (9* > 0). Tunnelling lifts their degeneracy, 
leading to energy levels _E± = Eq ± A. However, at degenerate bias the current vanishes, 
forcing one to move slightly away from this point. In order to find the levels for |/ x | <C 1 
we expand Eq. (|42|) near its minima, retaining linear terms in / x , (5 and quadratic terms in 

r/l 

X, 6*. Define 9 J as the minima, shifted due to f x and /3: 

^^W^i^^; (43) 

that is, the upper (lower) sign refers to the right (left) well. The potential energy (}4*2"|) then 
reads: 



+ A r /'£ +£ r/ # (44) 

E.T 



where the operator correction 9 = 9 — 9 J , 



^ = - a -2^ ±27T/ -^ (45) 



2a 



2aV4 



a 



1 V 4 2 



(47) 



Combining ()44j) and (|%2*j) in (jlHjl we obtain quadratic quantum mechanical Hamiltonian for 
the flux qubit in the left and right well near the degeneracy point: 



4(2ft + l) p . 2 , ^J. 2 



ft 



2/3 



+ 



2£ c ri 2 + EjB r 



+ [2E c n 2 x + EjA r ' l x 2 ] 
+ 8E c n e n„ + ^C" V 



(48) 



where 



(jr/l _ j_ 
~ 2a 



-l±27r/. 



2ft 2 - 1 



V4 



+ 



v/4 



ft z 



ft 



2ft 1 



2ft - 



4ft 2 



a 



(49) 



The first term in square brackets in (J48|) is the Hamiltonian of LC oscillator of the flux qubit, 
which is slightly modified by the last term in (|48jl . the next two terms in square brackets are 
oscillator Hamiltonians for the flux qubit variables, x an d 9, respectively. The interaction 
of the 9 degree of freedom with the qubit LC circuit is given by next-to-last term in ()48|) . 

Assuming the frequency (LC) -1 / 2 of the qubit LC circuit is much higher than the junc- 
tions frequencies Ej/H, Ec/h we neglect the interaction of the qubit variables, 9 and x with 
the qubit LC oscillator. This is equivalent to the averaging of Hamiltonian (J48j) over the 
ground state of the LC Hamiltonian. Therefore, for the qubit Hamiltonian we obtain: 



H qb = (H*' 1 ) = EjUf + + <^ 2 > + [2E c n\ + EjA^x 2 ] + 2E c n 2 e + E 3 W 



(50) 



where 



8E c Ej (2ft + 1) 







a 



1/2 



8(3 E c (2ft + 1) 



E, 



ft 



1/2 



Next we confine ourself only to the ground state of ()50j) in either of the wells. 



T r/l 



-e + EjXjf + Ej— 



Qr/l 



ftjjj 



r/l 



where 



hujf = Ej 



ag 



1 T 2k f. 



2a z 



x 



2V4a< 



(3- 



+ 



2a z 



r/l 



4ft 



(51) 



(52) 
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e V «5 \ 2(4a 2 -l) 3/2 M 4« 2 -lV 4 2 11 

The ground state wave functions in left (right) well are as follows: 



The tunneling between two wells lifts degeneracy yielding the well known result for eigenen- 
ergies E± = (e l —e T )/2 ± a/ (e l —e T ) 2 /A + A 2 , which was given above in Eq. (pTTj) . The Eqs. 
flEIJ), ffH^jl . allows us to calculate the numerical factor A in (|3*T|) : 



vr V 9 \V4tt 2 - 1 4a 2 - 1 

^ /a v/8(2a + 1) / 2a 2 - 1 + 4a 2 - 1 



g 8a V V4a 2 - 1 « 
The average current in eigenstates E± is calculated from (JT9*)l : 



S£± Lr A 2 (a,p) £. 



iq = ^ = ± jj /' - J ( 56 ) 

q <9$ x 2tt V£ 2 + A 2 

If the deviation from the degeneracy point fx = is significant (e >> A) then the current 

(JHHJl reduces to its local value in a particular well I q — > ±/c*A(a, g)/2ir = de r l l /d&x- 



A. The matrix elements of the current operator 

As is seen from Eqs. (J33J), ()35jh it is necessary to calculate the quantity _\dty + / d& x) • 
The calculation yields the following result: 

\ de A „ vr A 

+ / 3 77T^» (57) 



where 



d$ x / d<$> x 2A 2 r 4$ A £ 



n«) = : 7;\ + : : vi ( -7 + 7« 2 - 2 « 4 ) (58) 



(2a 2 - 1) a (2a 2 + 1) / 1 5 
4aV4a 2 - 1 + (4a 2 - l) 2 \ 4 + 2* 
The correction due to inductance (second term in r. h. s. of ()58Jl) is usually small, however, 
it is responsible for nonzero value of non diagonal matrix elements of the current operator 
in the flux basis. In the eigenstate basis Eq. ()33|) transforms to: 

de 1 



11 



-ea z + Aa x ) (59) 



where we neglect the correction due to inductance. 
In the flux basis Eq. transforms to: 

The Eq. (|5Uj) is the main result of our calculations. It shows that the current operator 
in the flux qubit is not diagonal in the flux basis as it is implicitly assumed in most of 
papers on the subject. The non diagonal term comes from the finite inductance of the qubit 
loop. Though for the usual qubit design this term is relatively small, nevertheless, it might 
give noticeable effects for larger values of j3 in the arrangements when two flux qubit are 
inductively coupled via a term MI1I2 in the Hamiltonian, where M is a mutual inductance 
between qubit's loops, Ji, J 2 are the current operators of the respective qubits. 

In conclusion, we show that the finite loop inductance of a flux qubit results in additional 
non diagonal term in the current operator in the flux basis. The result is important in the 
arrangements with magnetic coupling of two or more flux qubits. 
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